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Abstract. We consider periodic quantum Hamiltonians on the torus phase 
space (Harper-like Hamiltonians). We calculate the topological Chern index 
which characterizes each spectral band in the generic case. This calculation 
is made by a semi-classical approach with use of quasi-modes. As a result, 
the Chern index is equal to the homotopy of the path of these quasi-modes 
on phase space as the Floquet parameter d of the band is varied. It is quite 
interesting that the Chern indices, defined as topological quantum numbers, 
can be expressed from simple properties of the classical trajectories. 
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1. Introduction 

Topological quantum numbers are different from quantum numbers based on 
symmetry because they are insensitive to the imperfections of the systems in which 
they are observed. In some sense, topological properties are robust properties. 
Topological quantum numbers have become very important in recent years in con- 
densed matter physics, where measurements of voltage and electrical resistance can 
be conveniently expressed in terms of them. 

If you consider for example a simple model of non-interacting electrons moving 
in a two-dimensional bi-periodic potential V{x, y) subject to a strong uniform per- 
pendicular magnetic field Bz and a low electrical field E^ , the Hall conductivity axy 
of a given filled Landau electronic band turns out to be proportional to an integer 
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(1) (T^y = —C 

C is the Chern index of the band, describing the topology of its fiber bundle struc- 
ture iiiii- 

We investigate in this paper the value of C as a function of the potential V. In the 
limit of high magnetic field Bz, the above model is mapped onto the well-known 
Harper model: the potential V is considered as a perturbation of the cyclotron 
motion, and the averaging method of mechanics gives an effective Hamiltonian equal 
to the average of V on the cyclotron circles. We neglect the coupling between the 
Landau bands 0. For a high magnetic field, (hence for a small cyclotron radius,) 
this transformation gives an effective Hamiltonian Heff{q,p) ^ V{q,p), which is 
bi-periodic in position and momentum (the phase space is a 2D torus), and an 
effective Planck constant heff = hc/{eB). In this approximation, trajectories are 
the contours of H^ff ~ V. Furthermore, the expression of h^/f shows that the high 
magnetic field regime corresponds to the semi-classical limit. This model will be 
the starting point of our study in the next section. For the sake of simpHcity, we 
will denote h^ff by H in the sequel. 

We will restrict ourself to the case where 1/h = A^ is an integer (but this is not 



a strong restriction as discussed in section 3.1). The spectrum of these Harper- like 
models has then a finite band structure. To each band n = 1 ^ A^ is associated a 
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topological Chern index C„. In this paper, we calculate all these Chern indices in 
the semi-classical limit N ^ oo, for a generic Hamiltonian H{q,p) on the torus. 

The way we compute these topological indices is the following. First, we use the 
fact that the classical dynamics generated by H is integrable, to construct quasi- 
modes (the usual WKB construction) in section 4. An eigen-function |(/Jn(^i,^2) > 
of H in band n, depends on two quantal parameters related to the periodicity 
conditions of the wave function on the torus. A quasi-mode |V'n(^ii^2) > is a 
quantum state localized on a trajectory which is very close to |v?„(6'i,02) > in the 
semi-classical Hmit -^ ^ oo,|g||||] (except when tunnehng effect occurs). 

Because of tunneling effect, the eigen-function |i^n(^i,^2) > jumps from one 
quasi-mode to another one as (01,6*2) are varying and follows a closed path on 
the torus. By treating correctly this tunneling effects when it occurs, we obtain a 
good approximation of the bundle of eigen-functions |(^„(6'i, 62) > by quasi modes 
|V'n(6'i, 6*2) >, for every value of (6*1, 02)- This gives us the principal result of this 
paper: the Chern index C„ is equal to the homotopy number /„ of this path, on the 
torus, see eq.(p2|). We obtain the Chern indices for the whole range of the spectrum, 
and show that the non-contractible trajectories play a major role. 

This result is based on properties presented in the appendix ^ specific for the 
computation of the Chern indices. The analytical methods we use rest on properties 
of zeros of the Bargmann representation (divisors of a holomorphic section). This 
appendix is self-contained, and general results presented here could be useful for 
other purposes where topological Chern indices are involved. 

There is a constraint on the global spectrum |P0||^, that X)n ^"^ = +1- In 
section 5, we give an explanation of this in terms of classical trajectories. It is shown 
that the non-contractible trajectories on the torus are responsible of this fact. This 
non-trivial total Chern number is crucial for explanation of the Hall conductivity 
0. A natural question is then to determine which bands n have a non zero Chern 
index C„ in the generic case? We answer this question in section 0, where we explain 
how to compute the Chern indices of the whole spectrum, from the classical Reeb 
graph. 

Our approach is quite similar to Thouless et al.'s one |l| where they treated 
Hf{q,p) — cos{kq) + ecos{p). The difference is that they treated H^ within a 
perturbation approximation (with respect to e ^ 1), whereas we use a semi-classical 
approach. This enables us to treat any generic Hamiltonian H , and to describe new 
examples with non trivial C„, different from the Thouless et al.'s one. 

Our results extend previous work of the author |gj|,||l^, where the Chern indices 
were calculated for a restricted range of the spectrum where the classical energy 
level T,E = {{q,p)/H{q,p) = E} consists in two or three contractible trajectories. 
In these cases, it was shown that C„ is zero except for special configurations of the 
trajectories. 

We give a numerical illustration of these results in section 6. It is shown that 
our results, although derived in the semi-classical limit A^ — > 00, are also valid for 
quite small A^. This illustrates the robustness of the topological properties. This 
numerical example is also used throughout the paper to confirm our calculations. 

In section 8, we complete this work by giving a new (but equivalent) derivation 
of the interpretation of the Chern index C for conductivity, more intuitive than the 
usual one based on the Kubo formula [Q: we show here that C is equal to the mean 
velocity of a wave packet on the torus, when (0i, 6*2) move adiabatically. 
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2. Classical trajectories of a generic Hamiltonian 

A generic classical Hamiltonian on the torus is a generic C°° Morse function 
H{q,p) periodic in q £ M and p £ M, such that 

(2) H{q + l,p) = Hiq,p+l) = Hiq,p). 

In this section we investigate the topological properties of the trajectories of H, and 
introduce the Reeb graph ||lq|, which summarizes them. The aim of this paper is 
to point out the rules that give the Chern indices from the Reeb graph. 
A trajectory q{t),p{t) is a solution of the Hamilton equations of motion: 

dtq^ dpH{q,p), 
dtp = -dqH{q,p). 

The trajectories q{t),p{t) evolve on the plane IR^, but because of periodicity, 
they can be considered as trajectories on the torus Tqp, obtained by identifying the 
opposite sides of the cell [0, 1] x [0, 1]. 

2.1. Properties on trajectories. Because of conservation of energy: d{H (q{t) , p{t))) / dt 
0, each trajectory is included in an energy level S^; — {{q,p)/H{q,p) — E} of the 
Morse function H. Moreover, the velocity on such a trajectory is never zero, except 
when there is a critical point of H where dH — 0. 

We can easily obtain properties on the topology of J2e- 

1. For non critical value of E, 'Ee is a closed one dimensional curve, and each 
trajectory can be identified with a connected component of S^;. Therefore, 
each trajectory is a periodic orbit, an oriented closed curve on the torus Tqp. 
Its topology (homology) is characterized by two integers {ni,n2)eZ^ which are 
the degrees in the q and p directions of q{t),p{t) for i = ^ T where T is 
the period. (ni,n2) are relatively primes, except if ni =0 (or 712 = 0) in 
which case n2 = --l,0or+l. A contractible trajectory has homology type 
m = n2 — and is clockwise or anti-clockwise. 

2. A generic Hamiltonian has always non contractible closed trajectories with 
type (ni,n2) 7^ (0,0). In this case, every other closed trajectory is one of the 
three types (ni, 712), (— "-i, ~Ti2) or (0, 0), and each of these three types exists. 
(Example: consider a perturbation of H{q,p) = cos{n2q — nip)). 

3. For non critical value of E, each trajectory belongs locally to a family param- 
eterized by the energy E such that every trajectory has the same type. (The 
critical trajectories are bifurcations between these famihes). 

For a given generic Hamiltonian, by a good choice of the lattice generators on the 
plane {q,p), it is possible to deal with trajectories of types (0, 0) (0, +1) and (0,-1) 
only. We will adopt this choice in the sequel of this paper. 

Now concerning the critical points and the critical trajectories : 

1. There are three sorts of critical points: a local minimum, a saddle point or a 
local maximum. From Euler formula ([Q p.29) we have ^{mins)+^{maxs) — 
^{saddles) = (The Euler characteristic of the torus is 0). 

2.2. Example. We will give now an example in order to illustrate these statements. 
This example will be used throughout the paper. 

Consider the following Hamiltonian: 

(3) H{q,p)^Ho{q,p)+H,iq,p), 

Ho{q,p) = cos(27rg) + 0.1 cos(27rp), 

Hi{q,p) = P [exp (-100(9 - qo? - 10(P - P^f)] , 
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where P is a functional operator which sums the function in each cell and makes it 
periodic: 

The trajectories generated by H are shown on the figure ffl, for go = 0.45, po = 
0.45. There are one minimum F, two maxima A,D, and three saddle points B,C,E. 




0.4 -0.2 -0 0.2 0.4 0.6 0.8 1 1.2 



Figure 1. Trajectories and fixed points of Hamiltonian (|^). 

The figure || gives a schematic view of the families of trajectories, with their 
bifurcations, and homology types. It is called a "Reeb graph" TZ ||l^. By definition, 
it is the topological quotient space TZ — Tqp/T where T — {Y,e}e is the set of energy 
levels. The topology of this graph gives the relative location of the trajectories (the 
q variable here is only to indicate that the dashed fine can be obtained by a section 
at constant p in figure El) . 



(0,0) 
B 



(0,0) 



(0,-1) . 



(0,+ l) • 



(0,0) 



Figure 2. Reeb graph of H: a schematic view of the Trajectories 
types and bifurcations of Hamiltonian (|^). Each point of the dashed 
line corresponds to a non contractible trajectory. They form a one- 
dimensional family. Each point of a solid line corresponds to a 
contractible trajectory. 

The Reeb graph can be done for any given generic Hamiltonian. This figure will 
always have a similar apparency, with one transverse (dashed) line (from g = 
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to q ~ 1) corresponding to the non contractible trajectories, and branches from it 
(soHd lines), corresponding to contractible trajectories. Critical energies are at the 
extremities of these branches. The final value of the Chern indices will be read more 
or less directly from this graph in section 7. 

3. The Quantum Hamiltonian 

In this section we introduce well known material that will be useful for the sequel. 

3.1. The Hilbert space Tip of the plane . The Hilbert space corresponding to 
the classical motion in the plane phase space T* (IR) ^ M^ is 

Hp = L^{M), 

(the index P will remind us that we deal with the plane phase-space). 

We choose a symmetric quantization procedure: to the classical Hamiltonian 
H{q,p) we associate a self-adjoint operator H in Hp by 

H = 2. o'^"i'"2 6xp {i2TTniq) exp {i2'Kn2p) + hermitian conjugate. 
Where c„i„2 are the Fourier components of the function H{q,p): 

(4) H{q,p)= ^ c„i,„2exp(i27mig)exp(i27rn2p). 

ni,7i2£Z^ 

Since if is a real valued function, the complex coefficients c„^„2 must satisfy: 

We denote by Tg [respectively Tp] the translation operator by one period: Tq 
translates hy Q = 1 a wave function ip{x) and Tp translates its Fourier transform 

iPip) hy P ^ I: 

(5) fglPix) = ^J{X~1) 

(6) fpil,{p) = iPip-l). 
We may rewrite equations (||) and @ as: 

Tq = e-Kp{-ip/h), Tp = exp{iq/h). 
Quantum mechanically speaking, the periodicity Eq. (g) reads: 

(7) [H,fQ] = [H,fp]^0. 
To continue, we now have to assume that 

(8) [fQ,fp]=0. 
It is easy to prove that: 

TQfp^e-'/^fpfQ, 
hence (^) is fulfilled if and only if there exists an integer N such that: 

(9) N = ^elN*. 

with h = 2Trh. This hypothesis (||) can be regarded as a geometric quantization 
condition, which states that there is an integer number of Planck cells in the phase 
space. The semi-classical limit h -^ corresponds then to the limit N -^ +00. 

We now assume that hypothesis (||) is fulfilled. 

Note that hypothesis is not so restrictive, because any value of h can be well 
approximated by a rational number h — m/N. It is then easy to check that 
[TmQ,Tp] = 0. In this case we just have then to consider the torus phase space 
[0, mQ[x |0, P[ instead of [0, Q[x |0, P[. 
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3.2. The Hilbert space T-Lt{6i, O2) for the torus. According to the commutation 
relations (|^) and (g), the Hilbert space L?{1R) may be decomposed as a direct sum 
of the eigen-spaces of the translation operators Tq and Tp: 

deid92 



(10) np = L^{IR) = f [ht {91,02) 
'Ht{0i,02) = \\ip> such that 



(2^)2 










'{g 




- exp{i9- 
= exp(i0; 
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with {9i,92) € [0, 27r[2 related to the periodicity of the wave function under trans- 
lations by an elementary cell. The space of the parameters (^1,^2) has also the 
topology of a torus, and will be denoted by Tg. 

The space T-Ct{9i, ^2) is not a subspace of L^{1R), the space of physical states, it 
is a space of distributions in the x representation. 

It is well known that Ht{9i,92) is finite dimensional ||^. To see that, let IV' >€ 
Ht{9i,92). The Fourier transform of ip{x) is 02-Floquet-periodic of period P, so 
-0(2;) is discrete, it is a sum of Dirac distributions supported at points distant from 
h ~ 1/N from each other. Moreover, ipix) is ^i-periodic, hence ip is characterized 
by the N coefficients at the N Dirac distributions supporting points in the interval 
gG[0,l[. 

Explicitly, a basis of T-Ct{9) is given by N distribution denoted \j, 9 >: 



\j, 9 >= ipj/{x) = ^= ^ cxp (-i7ii6'i) 6{x - Qj - III), j = 1,... ,. 



1 / 92 

with9, = ^(^J + ^ 

Eventually we get: 

dimff Ht{9i,92)^N. 

Because of EqjR) the Hamiltonian H is block-diagonal with respect to the de- 
composition Eq.([lO|). The operator H acts on Ht(^i,^2) as a A^ x iV hermitian 
matrix, its spectrum <j{9i, 6*2) is made of N eigenvalues and TV corresponding eigen- 
functions: 

(11) H\ip,,{9,,92)>^E.,,{9,,92)\cpn{0i,O2)>, n = l,...,iV. 

For a given level n, assuming that En{9i, 02) is never degenerate V6', the energy 
level En{9i,92) forms a band as (6*1, ^2) G Tg are varying, and the eigenvectors 
\i^n{9i,92) > form a 2D surface in the quantum states space. But for a fixed 
(6*1, 6*2), and any A € ff, A|V'n(6'i, ^2) > is also an eigenvector. So for a fixed level 
n, the family of eigenvectors form a complex-line-bundle (of fiber = (T 9 A ) in 
the projective space of the bundle Ht -^ Tg. The topology of this line bundle is 
characterized by an integer C„ G Z, called the Chern index (||l5|, chap 1, p . 139). 

The definition of this topological Chern index is presented in section |5.lj . 

3.3. The Hilbert space 7ic(^2) for the cylinder. In the sequel it will be useful 
to consider an intermediate decomposition of Hp in terms of states periodic in 
momentum only: 

Hp = L^{1R) = Jhc{92)^ 

Uc{92) = in > such that fp|V' >= exp(i6l2)|V' >} • 

Each state in He (^2) is a sum of Dirac distributions supported at points g„ — 
jf{n+^) with n& Z . The weights belongs to P{Z). 
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Conversely: 



nc{02)= fnT{0ue2)dei. 



For each 02, the Hilbert space 7ic(^2)is associated with the cyHnder phase space 
denoted Cqp, where {q,p) € M^ is identified with {q,p+ 1). 

From a state \tpp >e Tip = L'^{IR) of the plane, we can construct a state 
\'4>c >G Tic{d2) of the cylinder by the operation: 



-foo 



(12) |V;c(e2) >= P.J^P >- Y. exp(-m02)T^|V'P > 



n— — oo 



From a state \il)c >^ 'Hc{^2) of the cylinder, we can construct a state |V't(^i, ^2) >^ 
?^t(^i, ^2) of the torus by the operation: 



+00 



(13) \M0i,O2)>^PeMci02)>^ J2 exp{-ine,)T^\^Pci02) > 



4. The quasi modes 

4.1. Construction of a quasi-mode. The W.K.B. method allows the construc- 
tion of a quasi-modes for integrable dynamics. See ||§], Lazutkin [g| p. 235. 
First, we recall the definition and general properties of quasi-modes. 

Definition 1. If F is a given closed trajectory with energy E, a quasi-mode with 
error e, is a quantum state \'ip > localized near F, which satisfies: 

with h = 1/iV going to zero, with {E — E) = o{h). The trajectory F is called the 
support of \ip > and noted 

F = SuppdV^ >). 

If we want to improve this error e (and have an error exponentially small with 
respect to h), we have to take into account tunneling effect between F, and F' the 
"closest" trajectory F' which has the same energy E. (|Q,||l^,||l^). We will come 
back on this point. 

Proposition 1. ||s|| 

The interval [E — e, E + e] contains at least one eigenvalue of H. 
Let Aq, = [E — a,E + a]. If Aq, contains only one eigen-value E* , with eigen- 
vector \ip >, then (for normalized vectors) 



(14) 3p e i?. 



\ip > -e'l'\i, > 



e 
< — 



consequently, < f\^ >^ as soon as e < a\f2 (for our purpose, we will say in that 
case, that \ip > is a "good approximation" of \ip >. This condition will be sufficient 
to construct a bundle of quasi-modes with the same topology, same Chern index). 

Let us now give properties of quasi-modes specific to our problem. The explicit 
construction of the quasi-modes can be found in |M|2Q|pi|. 



FREDERIC FAURE 



P * 



Figure 3. Quasi- mode localized on a trajectory F of type (0,0) 
in cell (0,0). 

Contractible trajectories of type (0,0). For a non critical contractible trajectory F 
of type (0,0), we first construct a quasi-mode on the plane, \ipp >G Tip. From the 
usual W.K.B construction, this quasi mode is localized on a single image of F, for 
example in the cell (0, 0). See figure ^. 

The energy E of this quasi-mode is given by the usual E.B.K. condition: 

S{E) = {k + l/2)h + o{h) , fce^, 

with S{E) the surface inclosed in the trajectory F. 

A quasi-mode |^c(^2) > on the cylinder can be deduced by eq.(p^. A quasi- 
mode I ^T (^1,^2) > on the torus can be constructed by eq.(|l3|). Note that the 
energy of these quasi-modes E doesn't depend on {91,62). 

Non contractible trajectories of type (0,±1). If F is a non critical trajectory with 
type n = (0,±1), F is closed on the cylinder Cqp and a quasi-mode |V^c(^2) >€ 
'Hc{62) can be constructed for every 6*2. 

|V'c(^2) > is localized on F, for example in cell 0, see figure 0. 




Figure 4. Quasi-mode locaHzed on a trajectory F of type (0, +1), 
in cell (0) on the cylinder. 



The energy E of this quasi-mode is given by the condition: 
(15) S{E) = {k-02/27r)h + o{h), k e Z, 

where S{E) is the surface included in the unit cell, on the right of the trajectory 

F (here oriented by increasing p). We justify this condition in the next paragraph. 

Note that the energy £'(^2) as well as the support F(02) depend continuously on 

A quasi mode |^t(^i,^2) > on the torus can be deduced from |V'c(^2) >by 
eq.m 
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For contractible and non contractible trajectories, an other quasi-mode with the 
same energy can be constructed in the cell n by 

(16) \i'cA(^2)>=T^\^c{02)>. 

From the W.K.B. construction, the error e of these quasi-modes is of order o{h°°). 
The error e is due to tunneling interaction with "the nearest quasi-mode of nearby 
energy" (see the next section). 

Critical trajectories. If Eq is a critical energy, quasi modes of energy near Eq (in the 
interval [Eq ± 0{h)]) and locaHzed near the critical trajectory, can be constructed. 

See iiiiiii. 

Justification of quantization rule eq\lq. When constructing a WKB quasi-mode on 
a closed trajectory, the quantization rule comes from the phase matching. The 
phase of the Fourier transform of the quasi-mode evolves like (— / qdp) along the 
trajectory q{t),p{t). For a trajectory T of type (0, ±1), the phase accumulated over 
a period is then ip = S/h = — Jp qdp/h, where S is the surface on the right of 
the trajectory F (oriented by increasing value of p). The periodicity condition ( |lo| ) 
requires that (p = O2 [2tt] . This gives the result ( |lq ) . 

4.2. The semi-classical spectrum without tunneling corrections. We now 

investigate the global properties of the quasi-modes, with respect to 02- 

Consider fixed values of {0i,92) G Tg. A W.K.B. construction as previously 
described, gives us TV quasi-modes \ipT{6i,d2) >n with energies En{di,02) over the 
whole energy range, with n ^ 1 ^ N. This is a semi-classical approximation 
(Tsc{di,02) of the spectrum (t{9i,92) defined by eq . (|l l[) . In fact the construction 
was made on Tic (02) and does not depend on 0i, so we will note it cTsc{02)- (For 
contractible trajectories, it does not even depend on 6*2). 

From the above construction, each quasi-mode |'0c(^2) > depends continuously 
on 02- So each quasi-mode belongs to a function ^ : ^2 ^ | "00(^2) >€ 'He (02) 
which is periodic with respect to 02 (up to a phase), with a period 02 — 2Trm, 
with m & N. (Otherwise, if the function ^ were not periodic, the spectrum asc{02) 
would be infinite for a fixed value of 6*2)- The corresponding energy £^(6*2) and the 
support r(02) are therefore also periodic with the same periods. 

If the trajectory F is contractible (type (0, 0)), the energy E and the support T 
are fixed. They do not depend on 02, and the period is 82 — 2tt. 

For a non contractible and non critical trajectory (type (0, ±1)), eq.(p^ tells that 
the surface 02 -^ S{02) on the right of F is a strictly decreasing function because 
dS/d02 = —h/2'!r. Then 02 -^ r(^2) moves from the left to the right on the set 
of non contractible trajectories on the cylinder Cqp. This function reaches all the 
non contractible trajectories, and jumps over a connected component of contractible 
trajectories, via a critical point . See figure ||. 

The surface S{92) is always a decreasing function, but the corresponding energy 
£^(6*2) decreases for trajectories of type (0, -1-1) and increases for type (0,-1), with 
the behavior shown on figure 0, if we plot E{02) for ^2 = ^ 27rm. Figure ra shows 
this function E(02) but folded in the interval 02 = ^ 2it. 

So we obtain: 

Proposition 2. 

The global spectrum 02 -^ (Jsc{02) is the union of one single function '^nc ■ 02 -^ 
\i'nc{02) >€ 'Hc{02), and many others {'^c)i=i,..,N, '■ 02 -^ \ipc{02) >G ^0(^2) such 
that 

The support Tnc of ^„c is on the non contractible trajectories. The period of ^„c 
is 02 = 2TT.Nnc, such that Nnc.h ~ Snc the surface occupied by the non- contractible 
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f (0,+l) (0,0) (0,-1) 




Figure 5. The dashed arrows indicate the motion of the support 
of a quasi mode, 02 -^ ^{(^2), over all non contractible trajectories, 
as 02 increases. 



E l\ 



(0,+ l) 



(0,-1) ,1 



non contractible trajectories 

critical trajectory 
contractible trajectories (0,0) 



e! 271 Q 

2 Oj 



Figure 6. The energy of the quasi- modes as a function of 92- 



trajectories. In phase space, this support 62 -^ r„c(02) is a closed cycle of homotopy 
1 on the cylinder Cqp. 

The support Tc of each function (^c)j is on a fixed contractible trajectory. Its 
period is 82 = 27r. The number Nc of these functions is such that Nc.h ^ Sc 
the surface occupied by the contractible trajectories. Note that N — Nc + Nnc and 

We now discuss the error of the spectrum (Jsc{02)- 

The error of each quasi- mode is e ^ o{h°°). There are crossings in the spec- 
trum asc(02) for discrete values ^2; (caused by the non contractible function when 



Enc{92) = En 



or Enc{d2) = Ec{92)), see figure pi For fixed 6*2, away from a 



small neighborhood of every crossing value 6*3, the spectrum asc{02) has N discrete 
eigenvalues and each eigen-value is isolated from the overs by an interval of length 
a^ e. We deduce from properties [^ that the corresponding quasi-modes |V'n(02) > 
(in fact their image in TCt by the mapping Pe2 ) , are a good approximation of the 
actual eigenvectors \ipni9 1,62) >, for any 9i. 

Let us remark that this semi-classical spectrum crsc(6'2) does not depend on 9i and 
is thus infinitely degenerate. This correspond to the invariance by the translation 
(|l6|). This is not the case in the actual spectrum a {9 1,92) because of tunneling 
effect. 

Figure II shows the numerical spectrum of our example H{q,p), eq.(g). From this, 
we deduce what the semi-classical spectrum should look like: fig (p^. Compare this 
with the semi-classical spectrum figure ^. 

<^ 30(^2) is not a good approximation of the actual spectrum near the crossing 
values ^2- To improve that, in the next section, we will take into account the 
tunneling effect at each crossing of energy. It causes a splitting of levels and we will 
recover N well defined bands of energy. 
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4.3. The semi-classical spectrum with tunneling corrections. Consider a 
neighborhood C/ of a value 6'2, where there is a crossing in the spectrum asc{d2) 
between two energies: £'i(f?2) = -E'2(^2)- See figure |6[ As said before, this crossing 
is caused two non contractible trajectories when £'„c(^2) = -E'„c(^2) '^^ by one non 
contractible trajectory and one contractible one: £'nc(^2) = Ec. The crossing due 
to two contractible trajectories Ed — Ec2 is not considered here, because it does 
not occur for a generic Hamiltonian H. See [|ll|||l§l for results treating this non 
generic case. From eq.(lq), these energies correspond to two families of quasi- modes: 
iV'lc.m >, IV'2c,n2 >G 'Hc{02) where ni G Z , n2 € Z . We suppose that these 
quasi-modes are localized respectively on two trajectories Fi and r2, in cells ni and 

To improve the error e of the quasi-mode say {il'lcni >, we have to take into 
account the tunneling interactions with other quasi-modes which have nearby en- 
ergy |'0Jc,m >with j = 1,2 and m G Z . See figure 0- This tunneling interaction 
comes from non vanishing terms Aj^m =< V'lc.nil-f^lV'Jc.m >• We will treat the 
tunnehng effect at the leading order, by keeping only the dominant term. Consider 
first the terms ^2,n2 with j = 2, and n2 € Z. We keep the n2 term with greatest 
modulus: |^2,n2l- This term is o{h°°) and describe the tunnehng interaction be- 
tween \iplc,ni >and the "nearest" quasi-mode \^2c,n2 >• See figure 0. It is clear 
that the nearest quasi-mode is n2 = ni or n2 — ni zL 1, because other values of 
n2 are in more distant cehs. The other terms are o[h°°) with respect to leading 
order term. Consider now the terms Ai^m with j = 1, and m G Z. The dominant 
term is for m = rii: Ai^m =< V'lc.jii |-H^|'01c,ni >= Ei{92). This is the energy of 
the quasi-mode. The others terms, which describe tunnehng interaction between 
iV'lc.m > and its images |V'lc,m > in other cehs, are o{h°°) with respect to 1^2,712! 
and so negligible. 



A r, 



l> 



\y 




m > 

nj-l 



"1 



m > 
"1 



m > 



Figure 7. Tunneling effect between quasi-mode iV'lni > on tra- 
jectory Fi, in cell rii, and quasi-modes |V'2„2 >• We keep only the 
dominant tunneling interaction, here due to |V'2„i >. 



We have obtained that the dominant correction to the quasi-modes \^lc,ni > is 
described by the 2x2 tunneling interaction matrix: 

^1(^ 2) A2^nA02) 

1,1 \ ^2,712(^2) £^2(^2) 



A{e2) = {<iAc,m\H\i>3c,, 



> 



where |V'2c,n2 > is the "nearest" quasi-mode to \iplc,ni >, ^2 = ni or n2 — ni± 1, 
and |A2,„2l is o{h°°). 

By diagonalizing the matrix A, we obtain two quasi-modes |i/)c-i-(^2) > and 
\ipc-{(^2) >G 'Hc{02) expressed as a linear superposition of IV'lc.ni >and |'(/;2c.ri2 >) 
and two energies £+{62) > £^-(6*2) with a gap \E+i92) ~ £-(^2)! > |-42,„2l • But 
the error of the quasi-modes is o{h°°) with respect to this gap. From property (|l|), 
we deduce that these quasi-modes |'0c±(^2) > (their image in TYt actually) are a 
good approximation of the eigenvectors |(^i(0i,02) >, for any 9i. 
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Now if we match these results with the outside of the neighborhood L'^ of 6*2 , we 
can obtain a continuous dependence of |^c(^2) >with respect to 62 £ M. 

Proposition 3. We have obtained quasi-modes |V'c(^2) > on the cylinder, which 



are a good approximation of the spectrum eq.(ll) for any 6*1, 6*2 in the sense that 



(17) < ipn{0i, 82)1^01,02) >^ for every (0i,02). 

We note this approximate spectrum by o'sct(^2) (the energies do not depend on 

There is no more crossing in the energy levels, so we have N well defined (semi- 
classical) bands and corresponding quasi-modes |V'c.ri(^2) >, n = l, ..,N. 

At each avoided crossing, we have obtained that the support of a quasi-mode 
jumps from the trajectory Fi in cell ni (resp. T2 in n2) to the closest trajectory T2 
in cell 7^2 — rii or n2 — ni zL 1 (respect. Fi in ni ) as in figure |j. 

More globally, for each band n, we derived an function for the support of the 
quasi-mode \ipc.ni02) >on the cylinder Cqpt 

(18) Sen-. 02 e i? -^ Supp {HCAO2) >) C Cqp. 

The image ST,n of Sc.ii on the torus is periodic for 612 = ^ 27r, and Sc,n 
is a lifting of ST,n on the cylinder. With Sc,n is associated a homotopic number 
In € Z. Precisely, In is given by: 

(19) T^" [Supp (I Vc.n(O) >)] = [Supp (|Va„(27r) >)] . 

/„ is simply the result of the jumps of the quasi-mode from cell rii to cell ni + /„ 
as 62 = 0^ 27r. 

5. Chern indices 

5.1. Definition . The eigen-functions [ipn{0i,62) > of the band n, eq.(|ll[), form a 
complex Hne fiber bundle over the torus {91,62) € Tg. On a contractible subset of 
U C Te we can choose the states [ipn{di,d2) > such that they are normalized and 
that they form a (continuous) section over U. Note that this is not possible over 
the whole torus Tg except if the bundle is trivial. 

The global topology of this complex line bundle is characterized by its Chern 
index (^ p. 139, ^). 

Because of the natural Hilbert scalar product on L'^{1R) — J J Ht {0i,92) dOi d92 , 
which induces the Berry (or Chern) connection ([g6|, [||), this topological number 
is explicitly given by the integral of the Berry (or Chern) curvature (|p]|,|[l^] p. 141): 

(20) C'n = -- / {<de^Lpn\de2^n> - <d0.,ipn\de^>^n>)d9id92, 

obtained by sum over local open subsets Ui which covered Tg, and sections chosen 
in each of them. 

This expression has been used intensively for our numerical calculations. 

Moreover, it can be shown (see e.g. [|o|,|g|) that: 

N 

(21) E ^" = 1- 

n=l 

There is an alternative expression for C„ more suitable for our analytical calcu- 
lations, given in eq. (p7|), (see [15| p. 141). It is based on the motion of the zeros of 
the Bargmann representation of the states [^n{9i,92) >■ 
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5.2. Chern indices of the semi-classical bands. Our aim is to express the 
index C„ from the description of the classical dynamics. 

Because of property ( p7| ) , we deduce from theorem (||) page g|,that the band of 
the semi-classical spectrum asc,t have the same topology (and same Chern indices) 
as the actual energy bands. 

Our work consists now in computing the index C„ of the semi-classical band n. 

From theorem m) page Ea, the result is simply that: 



(22) 



^n — -'ni 



where /„ defined by eq. ([Tg|), is the homotopic number characterizing the path 
followed by the support of the quasi-mode on the cyhnder, when 02 is varying from 
to 27r. 



5.3. Global analysis of the Chern indices. In this section we calculate the sum 
of Chern indices for consecutive bands, in order to recover the result eq.(|2l|) within 
our semi-classical approach. 

For each band n = 1 ^ A^, we have defined the cycle of the support ST,n on the 



torus Tqp, eq.(|l8[). The last result eq.(p2|), is that the Chern index C„ — I{Sn) is 
the homotopy number of this cycle (in the q direction) . (We will now drop the T 
suffix in Sn)- 

We decide to define the sum of two or more consecutive cycles Sn + Sn+i + 
... + Sn+a by removing the jumping due to tunneling at the crossings between two 
consecutive cycles Si and S'i+i. This is illustrated on figure @. The homotopy is 
then I{Sn + ■■■ + Sn+a) = HSn) + ■■■ + I{'^n+a), because the removing is a local 
operation. 



A) 





=> 




Figure 8. The tunneling effect generates a splitting between en- 
ergy levels (A), and a jump of quasi- modes between trajectories 
(B). In order to calculate the sum X)n=i '-'n, we have to remove 
this splitting, as well as this jump. 



^N 



J2n=i ^n is then the support of the semi-classical spectrum o-gc without tunneling, 
and from proposition (||) , it is composed of one cycle Tnc of homotopy 1 and many 
constant cycles Fci of homotopy 0. This gives: 



N N / ^ \ 

n— 1 n— 1 \n— 1 / 
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We recover eq. (pi|). The total Chern index (+1) is interpreted here as the motion 
of the non contractible quasi-mode on the torus as 62 is varying (in accordance with 

th.r 



page 



25). 



6. Numerical illustration 



We will illustrate the above results on the Hamiltonian eq.da). The numerical 



calculations of the Chern indices have been done with the curvature formula eq.(|20|) 
For TV = 11 levels, the Chern indices are: 



Ci^4 C5 Cg C-j Cs^ii 


+1 -1 +1 



Our analytical results have been obtained in the N ^ 00 limit. Although, A'^ = 11 
is not very high, the numerical results which follow can well be interpreted in the 
semi-classical description of this paper. 

The minimum value of N required to get correct semi-classical estimates, can be 
evaluated in the following manner: the mesh provided by Planck cells of magnitude 
1/N has to be about so fine that all phase space structures in a plot like figure |l| 
can be resolved. 

Figure^, shows the energy levels En{0i,62) for n=l— >ll,asa function of 62. 
The dependence on 9i gives a width of the levels, but is very weak and not visible 
on the figure. 
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Figure 9. Energy levels E„(^i, 02) for rt = 1 ^ 11, of Hamilton- 
ian eq.(Q), as a function of 02- There is no degeneracy, but the small 
splittings are not visible on the figure. Points A,B,C, correspond 
to the Husimi representations of figure |l^. 

In this energy spectrum, one can clearly distinguish two categories of energy 



levels. These two categories are reproduced on figure 10, where we have artifi 



daily dropped the splittings. This figure corresponds to the semi-classical spectrum 
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There are four energy levels (in solid line) which do not depend strongly on 02 
and have periodicity 2tt. In the light of prop. (||), these energy levels correspond to 
four quasi-modes |V'c(^2) > localized on contractible trajectories. There is also one 
continuous energy level -E„c(^2) (in dashed line) which depends on 62 and have pe- 
riodicity 27rm with m = 7. This energy level corresponds to quasi-modes \ipnci02) > 
localized on non-contractible trajectories. 



Figure 10. The semi-classical spectrum (Tsc{02) obtained from 
fig.(||). The four solid lines correspond to quasi-modes |'!/'c(^2) > 
locaHzed on contractible trajectories. The dashed line corresponds 
to quasi-modes |?/'nc(^2) > localized on non-contractible trajecto- 
ries. 

Figure O, shows the Husimi representation (see appendix (Kh) of eigenstates 
|</'n(^i,6'2) > for the level n = 6, with 6'i = , and three different values of ^2- 
These three eigen-states correspond to the points A,B,C on figures (^ and (p^. 
Looking at the classical trajectories on fig(yj), one can clearly associate these quasi- 
modes respectively with the trajectories: 

(A): with a contractible trajectory (point D on fig|l|), 

(B): a non-contractible trajectory of type (0, -1-1), 

(C): a non-contractible trajectory of type (0,-1). 



(B) 



(A) 



TT 



(B) 



! 



(C) 






(C) 



(A) ■ 



Figure 11. The dark density are the Husimi representation of 
eigenstates \ipn{0i,02) > for the band n — 6 (see figure ||), with 
6*1 = , and three different values of 02- (A): 6*2 = 0; (B): 
62 = 2,5, (C): 6*2 — 4,5. The dashed lines are the trajectories 
where the quasi-modes will jump at the next crossing. 



As 02 varies from to 2tt, the quantum state \'Pe{0i,02) > of level n = 6 jumps. 
From figure and figure M it is evident that there are three (avoided) level crossings, 
where the quasi-mode jumps: 
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Now from figure |lj, we can see that for the first crossing from (A) to (B), the 
quasi-mode stays in the same cell (Ani — 0). For the second crossing from (B) 
to (C) the quasi-mode changes by Ani = — 1 cell. For the third crossing from 
(C) to (A), the quasi-mode stays in the same cell (Ani = 0). This gives a total 
change of Ani — —1 cell and a homotopy /g — —1 for this sequence of quasi-modes. 
Accordingly, we have C^ — Iq = —\. Our result eq.(p2) is well verified here. 



7. The Chern indices of the spectrum in a generic situation 

From the above results, we give here the precise values of Chern indices for a 
simple and generic Hamiltonian (i.e. stable under perturbations) corresponding to 
the numerical example eq.(|^) and explain how to read them from the Reeb graph 

(fig- i- 

To simplify the discussion, we suppose in this section that we deal with a Hamil- 
tonian H such that the curve representing the non contractible trajectories in the 
Reeb graph, has only one maximum and minimum, like the dashed curve on fig. 
g. We suppose moreover that in the energy range of this curve, there is only one 
family of contractible trajectories, like the solid line (C-D) on fig. g. 

In the range of energy of the non-contractible trajectories, each band of the 
spectrum is made with one of the three following sequences of quasi-modes, see 
figure (H): 

• The band is made with the two non contractible trajectories of type (0, ±1). 
This gives the sequence: S = (-1 — ). 

• The band has also a quasi-mode on the contractible trajectory (c). This gives 
two possible sequences: S = {c^ — ) or 5 = (c — h). 
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Figure 12. The three possible band structures (solid lines), made 
with a non contractible trajectory of type (0, ±1) and a contractible 
trajectory (c). The dashed lines are the upper and lower bands. 



To simplify the discussion, we exclude sequences as {c^c—') or (c — c-t-). (These 
sequences do not appear if the energy intervals between successive energy levels of 
contractible quasi-modes are large enough). 

It is clear from figure |lj that if the band n has the sequence Sn = [c^ — ), then 
the band n -H 1 has necessary the sequence Sn+\ = (c — h), and conversely. These 
sequences come always by pair. 

To illustrate this from our numerical example, we can see on figure (||) and 
figure dlQ), that the bands 3 to 9 have respectively the sequences: S'3 = (c^ — ), 

^4 = (C - +) ,^5 = ( + -) ,^6 - (C + -) ,57 = (C - -f ) ,^8 = (+-), ^9 - ( + -). 

Now we saw in prop. (S) that a non zero Chern index results from the jumping of 
neighboring quasi-modes on the cyHnder, via non contractible trajectories. We saw 
also that a given quasi-mode IV" > jumps by tunneling effect to the nearest quasi- 
mode of the same energy located on a non contractible trajectory. It is therefore 
important to determine in which direction (right or left) this jump occurs: we have 
to complete the Reeb graph (fig 0), by giving the direction of the shortest jump 
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between any two trajectory. For that purpose, it is sufficient to give the location of 
the separatrix between left-jumps and right-jumps. 

On figure |l^, we have reproduced the Reeb graph (the solid curve) , and we have 
drawn (dashed line and dotted line) the possible location of this separatrix. The 
effect of this separatrix (dashed line) is shown with the plot of the jumps from one 
particular quasi-mode 1-0 > to a non-contractible trajectory of type (—1) at two 
different energies, on both sides of the separatrix. The separatrix for the jump 
from/to a non contractible of type (-1-1) is in dotted line. 

El is the energy where this separatrix crosses the non-contractible trajectories 
family. That the energy Ei is the same for the dashed and the dotted line comes 
simply from the reciprocity property of the tunneling "geodesic". 

What matters in fact is only the intersection of this separatrix with the Reeb 
graph. The picture |l^ is qualitative here, but quantitative calculations of the tun- 
neling effect would precise the energy of these intersections. 




Figure 13. Example for the jumping of a quasi-mode |'0 > to 
the nearest non contractible trajectory (of type (0, —1)). The solid 
curve is the Reeb graph as in fig. (||). The dotted and dashed lines 
represent the separatrix for these attraction basins. 



Now, by combining figures |^ and |1J together with rule eq.(g2|), it is easy to 
compute the Chern index of the three kind of bands. This is done on figure |lj for 
the sequences S'5, 5*6, S's, and we infer the general result: 

• If El is in the range of a band of type (H — ) then its Chern index is C(+_) = 
-1-1, otherwise C(+_) = 0. 

• If the trajectory (c) is in the region R (of fig. |l^ then C(^c^ ) — —1 (respect. 

-1-1), if the band has energy greater than Ei (respect, lower). In the other 
regions, C(c+-) = 0. 

• If the trajectory (c) is in the region R then C(^c h) = +1 (respect. —1), if 

the band has energy greater than Ei (respect, lower). In the other regions, 
C(c-+) == 0. 

This allows us to interpret the Chern indices of our numerical example: the band 
^5 = (H — ) with C5 = -1-1 corresponds to the energy Ei. Above it, the two bands 
5'6,S'7 with Chern indices Cg = — 1, C7 = +1 form a pair, thanks to the presence of 
the contractible trajectory (c) in region R. 

As a general result, we conclude that for any choice of Hamiltonian H with a 
similar Reeb graph, there is a unique band in the middle of the spectrum, including 



FREDERIC FAURE 




C=-l 



Figure 14. From knowledge of the sequences S5, Sq, Ss, their ho- 
motopy and the Chern indices C5 = +1, Ce = — 1 and Cs — are 
graphically computed. C5 = +1 is associated with a band of energy 
El. Note that Ce = — 1 is non zero thanks to the quasi-mode (A) 
on the contractible trajectory (c) in region R. 

energy Ei, with Chern index +1. Above this energy (for Ei < E < E2), the Chern 
indices come by pairs (— 1,+1) , and under the energy Ei, the Chern indices come 
by pairs (+1, —1). 

Of course for another appearance of the Reeb graph, the rules could be different. 
In particular, the Chern indices can reach higher values if the non-contractible curve 
of the Reeb graph has more than one maximum. We do not know if there a simple 
and explicit expression of the Chern indices in the general case. 



8. Chern index and Hall conductivity 

It is well known |l| that the Chern index C of a band is related to the transverse 
Integer Quantum Hall conductivity axy by the relation eq.(|l|). The usual demon- 
stration of this formula starts from the Kubo linear response formula for axy, and 
expresses it as the curvature integral for C, eq.(po|). In this section we want to 
demonstrate it from an other point of view, by relating C with the motion of a 
wave packet on the plane Tqp. This result seems more intuitive. 

Consider a given energy band n of the spectrum eq.([Ll|). For each (^i, 62) there 
is one eigenvector |(/J„(^i, 6*2) > defined up to a multiplicative constant. Suppose 
that we have a C°° function (6li,6l2) ^ |V'(^i, 6*2) >G HTlfi*!, 6*2), where |V'(6'i, 6*2) > 
is proportional to |(^n(^i, ^2) >, eq.([ll]). Note that if C„ 7^ 0, this function must be 
zero at some values of {0i, 62). 

Consider now 



IV' >= 



deide2\'ii^{ei,B2) > 



Then \ip >^ Hp = L'^{M). We can say that |?/' > is a state "locaHzed" on the 
plane Tqp constructed from the band n, like Wannier states in Solid State Physics. 
(Conversely, given |$ >€ L^(_ZR), we can project it on the spectral band n and 
obtain such \ip >■) 

The state \^ > can be translated by {ni,n2) cells on the plane by: 



l^n„n. >=r5^r;=IV^>= 



^6*1^6126 



ini0i+in2O2 



moi,92)>. 
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Conversely: 

These simple relations can be generalized by the following result: 
Consider an other state constructed from the band n: 

(23) !(/)>= / [ deide2e'f^'''''^\ij{ei,e2) 



where / is a continuous and periodic function such that: 

/(0i+2^,02) = /(0i,02)+iVi27r 

/(01,02+2^) = /(01,02)+iV22^, 

{Ni , N2 e Z^ are topological integers which characterized the homotopy of the 
function /). 
Let us define: 

<ni >= ^ ni. |< ■tpm,n2\4' >l 

ni,n2 
< n-2 >= ^ "2- |< '0ni,n2l0 >l^ i 

"1,"2 

which give the mean position of the state \(j> >on the plane Tqp relatively to the 
localized states \'ipni,n2 >■ It is easy to prove that: 

< ni >= A^i 

< n2 >= N2, 

the mean position of \(j> >is then "quantized". 

(This is a simple property of Fourier series: if g{6) = J2n CnC™^ = e*-'^^^' is 
periodic and f{2n) — /(O) + N2tt then TV = X^n^l'-"!^- Another way to say it is: 
X^ri'^'I'^nP — < 9\pe\9 >— -^ /o f'd9 — N, with the current operator pe — id/dd.). 

We now give some consequences of this result on the "quantized mean position". 

Now let us first look at the time evolution of |^o,o >: 



\^(t) >= e-'^'/''|V'o,o >= / / d9,d92e-^''^'^''^^'/^\^{e,,92) > 



Comparing to eg. (|23|) , the phase function is the dynamical phase / (6*1, 02 ) = —E{9i,92)t/h 
of homotopy type A^i = A^2 = . (because by continuous homotopy t ^ 0, 
/ is mapped to 0). This means that during its evolution, the quantum state 
\tp{t) >spreads over the plane, but its mean position < ni >,< n2 > stays zero. 

Let us suppose now that 6*1, 612 are no more "good quantum numbers", and that 
there is a slow drift. In the model of bi-dimensional electrons of the introduction, 
this drift is the adiabatic motion due to the low electrical field E^. Consider for 
example: 

9i{t) = -ujit 

6*2(0 = 6*2(0). 

After one period T — 2tt/uji, the evolution gives 

e-'"^^^\iP{9i,92) >= exp {icJ)d{9i,92) + i(l)Bi92)) m9i,92) >, 

where (/'d(6'i, ^2) is the dynamical phase of homotopy type A^i = iV2 = and 
0b(^1j 6*2) is the Berry's phase of the path 9i{t) and doesn't depend on 9i henceforth. 
By homotopy deformation, </'b(6'2) = 27rC02 and is therefore of type A^i = 0, A'^2 = 
C This means that after one period T, the mean position of the quantum state 
\ip(t) >has been changed by i5 < ni >= 0, S < n2 >= C (in cells units). 
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This "quantized" velocity V2 — S < n2 > /T is responsible for the Integer Hall 
conductivity. Indeed, in the Harper Model [Q of non interacting bi-dimensional 
(x, y) electrons in a bi-periodic potential of period {X, Y), and a perpendicular high 
magnetic field Bz, with a low electrical field along E^, the adiabatic hypothesis 
gives: 

h 

We have obtained that the quantized velocity of one electron quantum state is 
Vy = {CY)/T. For a filled band, the electronic density is one electron per cell: p = 
1/{XY). The electronic density current of a filled band is then jy = peVy = j^C.E^ 
which leads to eq.(|l|). 

9. Conclusion 

By using quasi-modes, we have obtained a semi-classical description of band 
eigen-states of a generic quantum Hamiltonian on the torus phase-space. By using 
general results exposed in the appendices, relating the topological Chern index of 
a band, with the localization properties of the quantum states, we have been able 
to express the Chern index of each band in terms of the motion of the quasi-modes 
on the classical phase space. 

In section 7, we have shown how the Chern indices of the whole spectrum can 
be graphically computed from the classical Reeb graph. There is an interesting 
question which can be raised from our work: conversely, from a given sequence 
iCn)n of integers such that X)n ^» = +-'-; ^he question would be to determine 
the potential V{x,y) whose band spectrum possess this precise sequence of Chern 
indices. 

Currently one tries to observe experimental signatures of the Harper spectrum 
(Landau level substructures) in superlattices with periods of about 100 nm on GA- 
As-AlGaAs heterojunctions (|p7||,||2^). In these experiments, the Chern index gives 
the quantum Hall conductivity axy- In this paper we have shown simple corre- 
spondances between the Chern index and the classical trajectories of the electrons. 
This correspondances could be observed in experiments by measuring the varia- 
tions of the integer Hall conductivity, in relation with the potential V{x, y) of the 
supperlattice, which may be created by an external electrostatic grid. 
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Appendix A. Bargmann and Husimi representation 

In this appendix, we recall well known results on Bargmann and Husimi repre- 
sentation on the torus phase space (||l9|). 

We have seen previously that the space Tix is not a subspace of L'^{IR). For 
functions belonging to Ht, it will be more useful to introduce a phase space repre- 
sentation of the quantum states, called the Bargmann representation ( ||29|| ). 

Consider a quantum state |V' >G L^{M). In order to characterize the localization 
of IV' > in the phase space near the point (g,_p), we first construct a Gaussian wave 
packet \qp > (coherent state) defined in the a;-representation by: 

/ 1 \ ^^■^ i (x- qf 

The notation \qp > recalls that the coherent state is localized (in the semi-classical 
limit) at the point {q,p) of the phase space. 
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The Husimi distribution of a state |V' > is defined over the phase space by: 

and for ip £ L'^{]R), we have: 

(24) ||^(,)|2rf, = |||<,p|^>|2^. 

To characterize the functions of L^{T*1R) which are {q,p) representations of a 
state, it is more convenient to introduce a complex representation of the phase space 
z = -^{q + ip)- Another (proportional) expression of the coherent state is then: 

\z >= exp(za+)|0 >, 

with |0 > being the fundamental of the harmonic oscillator Hq = q'^ + jP, and 0+ 
being the associated creation operator. Indeed: 







1 ^ f.qp q^+p^-,, ^ 






The following 


anti- 


holomorphic function of z is called the Bargmann 

b^{z)^<z\^>. 


distribution 


of 


Clearly, we 


have 


h^{q,p)^\h^{zfe ^^, 







hence the zeroes of the function h^{q,p) are those of the holomorphic function 
6i/,(z), which are localized zeroes in the phase space. Moreover, ( |24| ) implies that 
ip e L'^{1R) if and only if 6^ G L'^^^ e"'^' /'') and h^ is anti-holomorphic. 

The same definitions can be applied for a state |V' >G Wt(^i,^2)- The corre- 
sponding Bargmann function is a theta-function [gO| and the Husimi distribution is 
bi-periodic in (q,p), hence is well defined on the torus Tqp. 

We need the following properties concerning the Bargmann representation b^ (z) =< 
zW: > of a quantum state 1-0 > on the torus (|z > is a coherent state on the torus 

If IV' >G 'Ht{0i,92), the Bargmann function b-^{z) is (a theta-function) anti- 
holomorphic with respect to z, and with N zeros Z = {zi',...,zn) in the cell 
[0,1] X [0,1]. These zeros are constrained by ([Q, this is "Abel theorem", and 
"Jacobi inversion theorem", see ||5| p. 235): 

(25) t^.^"i 

i=l 

Efl 

i=l 

with Zi — qi + ipi. Conversely, such a collection of N zeros define a state \ip >& 
'Ht{0i,92), unique up to a multiplicative constant. 

Note that since the N zeroes are constrained to have a fixed sum, we get the 
right dimension iV — 1 + 1 for the Hilbert space T-Lt{Oi, 6*2). 

Appendix B. Calculation of the Chern index in special cases 

For simplicity, we will note — {61,92) £ IB?. 

We will compute the Chern index of a line bundle in particular cases. We will 
always consider a line sub-bundle of TCriO) -^ Tg. Such a line bundle F is charac- 
terized by giving local section 9 ^ \ipii(^) >G 'Ht{9) over open sets Ui which cover 
Te. 
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From appendix ^ the line bundle F is also characterized by the function 9 — > 
Z{9) , where Z{9) is the set of N non-ordered zeroes of the Bargmann functions 
^V(e)(^)- 
B.l. Theorem on homotopy invariance. 

Theorem 1. . Consider two complex line bundles over Tg with Chern index C and 
C" defined by local sections \'ipi{9) >and IV'K^) >on open sets Ui. Suppose moreover 
that 

< MO)mo) >^o, y9, vi, 

then C :=:€' . 

Proof. We can suppose that the sections are normalized (< '!/'i(^)IV'j(^) >= !)■ On 
a set Ui, and for a given 6*, define: 

pm = \<uo)wm>\<^. 

and for fixed A e [0, 1]: 

\^,^^{9) >= \M0) > ■ < MO)mo) > +x{m9) > -\Me) >< M^Mm >)■ 

In the theorem, we suppose that Pi{9) > 0. Thus 

< ^Ii9)\cp,^^i9) >^ p^ + Xil - p^) > 0, 

and then \tpi,xi9) >^ 0. 

We now check that the sections \'fi^\{9) > define a complex line bundle F\. On 
an other set Uj, if |'0j(6') >= exp(ia) \tpii9) >and \i'j{9) >= exp(?/3) \^'^{9) > then 
\fj,\{(^) >— exp(j/?) \(pi^x{9) > is in the same line than \ipi^\{9) >. 

We have thus obtain a homotopic deformation F\, A S [0,1] between the line 
bundle Fq and Fi . They have thus the same Chern index C = C\ ~ C . D 



Theorem 2. More generally ( jlg^ ,p. 14-1) if we suppose that 

9^<^,{9m{9)>, 
has zeros in 9* with index i{9*) — ±1 ,then 

(26) C = C'+ Y^ i{9*). 

zeros 9* 

B.2. The Chern index from the zeroes of the Bargmann function. Let F 

be a line bundle as above. 

• Suppose that there exists some point zq S Tqp of phase space, such that 

(equivalently, this means that zq ^ Z{9), for all 9). Then we can select a 
vector \ip{9) >G H]sj{9) in each fiber such that arg (6^(g)(zo)) = 0. This gives 
a non- vanishing global section of the bundle. This is also a global frame, hence 
the bundle is trivial, and C — Q. 

• More generally, define 

7V(zo)-{0eTe / zoeZ(0)}, 
then 

(27) C= E (±1)' 

eeN(za) 

where the sign ±1 corresponds to the local orientation of the mapping Zi at 
9, where Z = {zi, . . . , zjv} and Zi(9) = zq. This can be deduced directly from 
eq.(||), with |-0' >= |zo >• 
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B.3. Example of a line bundle with a given Chern index C. We now con- 
struct explicit examples of such a bundle. 

• Suppose that N > 2. We define a line bundle Fq, by specifying the zeros of 
the Bargmann representation of a section, with the notations of eg. (|25|) : 

(28) 91 "fe' P'^^ 

92 = 0, p, = - 
q^ = 0, p, = 0ior3<i<N. 

The constraint eq.(pa) is easily verified. If we choose zq = (go + *Po) with 
Qo — Po — 0.5, Zq ^ Z{6) and from eq. (p7|), we deduce that the bundle Fq is 
trivial, with Chern index C = 0. 

• More generally, consider the bundle Fc defined by 

qi^il-C)^, pi=0 

q,^0, p,=0 ioY3<i<N. 

The constraint eq.g^is verified. We can apply eq. (P7|), to calculate the Chern 
index C . The set N{zq) is given by 

e,=TT[27r], 92 = ^[27t/C], 

and we obtain that the Chern index is C. 

B.4. Chern index for a moving coherent state. 

Theorem 3. Suppose N > 2. Consider the line bundle Fz defined by the local 
section \tpid) >= \ze >& Ti-TiQ) ^^'^ zg = qg + ipg, and 

nil ni2 \ ( 0i/27r 



(This means that the coherent state \zg > is moving over the cells as 9 is varying.) 
The Chern index of this bundle is then 



(30) C^N. 



nil ni2 
n2i n22 



n2i + ni2. 



Proof. Consider the fine bundle Fq defined by the zeroes eq. ( |2^ ) , with local section 
noted \(pi{9) > on each open set Ui C Tg. This bundle has C = Chern index. 
Consider the function 

f : e ^< ip^i9)\zg > . 

From eq.(p^), the Chern index C of the bundle Fz is 

c= Yl ^(^*)- 

^^zeroes of f 

The zeroes of / are given by: 

< LPi{6)\zg >— <^ Zg IS a zero of the section \ipi{9) > 

^h'='^' ,i=l,...,7V. 
{ Pe=Pz 

this gives 

nii^+ni2^ / |J ( / . . „ ., 

X , ll = \ n '^^ = \ 6i or = <^ for z = 3, . . . , N. 
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deduce from eq.(|) that 






C = 


"11 "12 - 1 
"21 "22 


+ 


"11 "12 
n21 — 1 "22 


+ (TV - 2) 


"11 

"21 


"12 
"22 


= 




N 


"11 "12 
"21 "22 


+ n21 + "12- 







D 

B.5. Bundle constructed from periodic motion of states on the plane. A 

natural question is : is it possible to generalize the theorem ||, for the periodic 
motion of arbitrary states in L^ (M) rather than only coherent states? 
There is a first result: 



Theorem 4. 

such that : 



Let S : 9 ^ \tp{d) >e L'^{]R) for e IB? he an arbitrary mapping 






(31) ye, \ij{0) >^ 

\ijiei+2TT,92)>-- 

|V(0i,e2 + 27r)>-.Q .p 

with ("11, "12, "21: "22) G ^^ which characterizes the periodic motion of the states 
\^p{0) > as 9 is varied. 

For a contractible open subset U ofTg, define 

\i^{9) >= Pe,PeA\^{9) >) ^'Ht[9). 

Then if N > S, and for a generic map S, \'4'{9) > gives a local section of a well 
defined line bundle over Tg, which can be noted also S. 

Proof. For every 9, the line bundle is defined by the vector \'4'{9) >. We have 
therefore to check that \ip{9) >^ , and that this line is periodic with respect to 9. 
Because the space Ht{9) is N dimensional, the condition \ip{9) >— is A^- 
dimensional. As soon as A^ > 2, and for 9 € [0, 27r]^the condition \^{9) >— 
cannot be satisfied generically. Now 

|V;(0i + 27r, 02) >= Pe,+2.(|V(^i + 27r) >) 
= Pe,Tl^''T';''\i^{9i) > 

= exp(mii6'i)exp(mi26'2)|V'(^) > • 

a2n, 92 + 627r) >are proportional and define the same line 

D 



So \^p{9) > and \ip{9i 
V(a,5)e Z^. 



Theorem 5. Suppose that N > 4. Let S : 9 ^ \%l){9) >& L'^{IR) be a generic 
mapping as defined in theorem Qa. Then the Chern index of the line bundle S is 
given by formula p3j. 

Proof We note 5*0 = S, and 5*1 = F^ the line bundle defined in theorem^. The proof 
consists in constructing a homotopic deformation Sx from 5*0 to 5*1 for A € [0, 1]. 
Because the Chern index C\ oi S\ is constant, we therefore conclude that the Chern 



index of S is given by formula (30). 
For A e [0, 1], and 9e R^, define 

\M0) >= MH9) > +(1 - A)l^(e) >e L\IR), 

where \z{9) > is a standard coherent state on the plane, z{9) — qg +ipg, with qg,pg 
given by eq. (|29|). We suppose moreover that the phases are chosen such that eq. 
(|l|) holds for \z{9) > and \i^{9) >. 

Then \ipx{9) >= if and only if A = 1/2 and \ip{9) >= -\z{9) >. This last 
situation is non generic, and if it occurs, we can just choose z{9) = {qg + go) + 
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i{p9 + Po) with arbitrary qo^Po to avoid this. So S\ : 6 —^ |^a(^) >G L^{1R) fulfills 
theorem |[ and defines a line bundle for every A. Because there are three parameters 
61,62, X we have to suppose now that A^ > 4 so that |V'(^) >¥" 0- □ 

B.6. Bundle constructed from a periodic motion on the cylinder. In sec- 
tion Q, we obtain the periodic motion of quasi-modes on the cylinder Cqp phase- 
space. A slightly different result than theorem || is then needed: 

Theorem 6. Suppose N > A, and let 

62e]R^\^ci02)>enc{62), 
be a generic continuous mapping such that 

V^2 G M, \iPc{62) >^ 0, 

with I <^ Z . 

For 6^ {61,62) define 

\^j{6) >= PeA^c{02) >(^nT{6). 

Then 6 -^ \''P{6) >€E Ti-rid) is a local section of a well defined line bundle over Tg, 
with Chern index I . 

Proof. The proof is similar to that of theorem ^. 

The bundle is well defined because, first \ip{6) >— needs N conditions, generi- 
cally not realized for arbitrary 6*1, 02 € [0, 27r], as soon as A^ > 2. 

Secondly, \^p{6) > and \i){6i+a2TT,62+b2Tr) >= Pe,T^^\ipc{02) >= exp{iIb6i)\tP{6) > 
are proportional and define therefore the same line. 

In order to calculate the Chern index, consider the mapping 

62eM^\zi62)>eL^iM), 

where \z{62) > is a standard coherent state on the plane, 2^(6*2) = qe + ipe, and 

qe = 1.62 

Pe = Po = constant, 

we can therefore construct \z{62)c >= ^6*2 1 -2(^2) ><= 'He (62) and \z{6i,62) >— 
Pe,\z{62)c>&nT{6). 

From theorem ||, 6* — > \z{6i,62) > define a line bundle with Chern index 7112 = /. 

Consider now for A e [0, 1] 

62 -^ \MG2) >- A|VXe2) > +(1 - \)\z{62)c >e nc{62), 

and 

Sx:6= {61,62) ^ \M0) >- Pe,\M02) >e Ht{6). 

Then for fixed A, S\ define a line bundle over Tg, because \ip\{6) >— nor \'ip\{6) >= 
are generic (cf discussion in proof of theorem ||), and there is also periodicity with 
respect to 6. 

The mapping A ^ S'a is a homotopic deformation from the line bundle 6 -^ 
I '0(0) > to the line bundle 6 -^ \z{6i,62) >■ Their Chern index is constant and 
equal to /. 
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